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ABSTRACT 

We study the dynamics of continuum driven winds from rotating stars, and develop an approximate analytical 
model. We then discuss the evolution of stellar angular momentum, and show that just above the Eddington 
limit, the winds are sufficiently concentrated towards the poles to spin up the star. A twin-lobe structure of 
the ejected nebula is seen to be a generic consequence of critical rotation. We find that if the pressure in such 
stars is sufficiently dominated by radiation, an equatorial ejection of mass will occur during eruptions. These 
results are then applied to ry-Carinae. We show that if it began its life with a high enough angular momentum, 
the present day wind could have driven the star towards critical rotation, if it is the dominant mode of mass 
loss. We find that the shape and size of the Homunculus nebula, as given by our model, agree with recent 
observations. Moreover, the contraction expected due to the sudden increase in luminosity at the onset of the 
Great Eruption explains the equatorial "skirt" as well. 

Keywords: stars: rotation — stars: winds, outflows — stars: individual (?/-Carinae) 



INTRODUCTION 

jy-Carinae epitomizes the class of Luminous Blue Variable 
(LBV) stars. During the 1840's Great Eruption, it shed sev- 
eral solar masses while radiating at several times the Edding- 
ton luminosity (Smith & Frew 2011). The Homunculus neb- 
ula which formed exhibits a clear bipolar structure with an 
equatorial "skirt" (e.g., Davidson & Humphreys 1997). 

Several ideas have been suggested to explain the broken 
spherical symmetry. One type of mechanisms relates the 
ejected nebula to the rotation of ?7-Carinae. A second type 
links the nebula to an interaction with the companion star. 

Frank et al. (1995) described the Homunculus as the re- 
sult of a fast nearly spherical wind interacting with a slow, 
previously accelerated, equatorial wind from a rotating star. 
Frank et al. (1998) describes an opposite model, where a fast 
aspherical wind interacts with a spherical one. 

Owockiefa/. (1998) and Maeder (1999) pointed out that 
von ZeipeFs theorem (von Zeipel 1924) has a significant ef- 
fect on the geometry of the wind. In particular, mass loss from 
a fast rotator is expected to be predominantly concentrated to- 
wards the poles, where the effective gravity and radiative field 
are large. This naturally explains the bipolar structure without 
requiring a two winds interaction, as previously suggested. 

Maeder & Desjacques (2001) applied the standard line 
driven winds model (Castor et al. 1975) to explain the form 
of the nebula. They have shown that if the polar temperature 
is ~ 2.5 x 10 4 K, the line driving bi-stability jump would give 
rise to a large increase in the equatorial mass loss, thus offer- 
ing an explanation to the origin of the "skirt". 

An alternative model for the overall structure of the Ho- 
munculus was developed by Smith & Townsend (2007). The 
nebula is described as explosively ejected material from a crit- 
ically rotating surface. An assumption that the initial velocity 
is proportional to the effective gravity gives the bipolar struc- 
ture. The skirt is obtained within some range of parameters, 
for which the orbits of matter originating at opposite low lati- 
tudes pass through the equatorial plane and collide. 

Recently, Soker (2004, 2005) advocated the idea that the 
twin lobe structure is due to an interaction with the compan- 
ion star. In particular, Soker argued that it is hard to reconcile 



the large bipolar mass loss of the Great Eruption with stan- 
dard models of stellar winds, for two reasons. First, one can- 
not obtain a high ratio between the wind momentum (M w Voo) 
and the photon momentum (L/c), as observed. Second, in or- 
der to explain the bipolar structure, the envelope is assumed 
to be near critical rotation. However, this is inconsistent with 
large mass loss, as it would slow down the envelope consid- 
erably. Soker also presented a model in which a slow wind 
from 77-Carinae is accreted onto the secondary, accelerating 
twin jets which form the Homunculus. As we shall see below, 
a model of super-Eddington (SED) continuum driven winds 
from a rotating star circumvents both problems. 

Shaviv (2000) has shown that the SED state can be ex- 
plained by a reduced effective opacity resulting from a 
"porous" atmosphere. In subsequent work, it was shown 
that such atmospheres give rise to optically thick winds 
with a relatively simple mass-loss luminosity relation, which 
also explains the behavior of classical novae (Shaviv 2001b; 
Owocki et al. 2004). The bipolar structure of the Homuncu- 
lus is a clear incentive to consider that its origin is a SED 
continuum driven wind from a rotating star. 

This note is organized as follows. In §1, we develop the 
theory behind continuum driven winds from rotating stars and 
the associated evolution of angular momentum. In §2, this 
theory is then applied to the case of 77-Carinae, explaining its 
high spin, and the formation of the Homunculus nebula and 
equatorial skirt. 

As a side note, we study in §A generic interactions be- 
tween the binary components, and show that spin-orbit cou- 
pling could not have been responsible for the high spin. In §B, 
we show that during the first few years of the Great Eruption, 
77-Carinae could not have been rigidly rotating. 

1. PART I -THEORY 

We begin in §1.1 by developing an approximate analytical 
model of continuum driven winds from rotating stars. In § 1 .2, 
we show that at low mass loss rates, such winds are capable of 
spinning up the star, even though angular momentum is lost 
to the wind. In §1.3, we study the consequences of eruptions 
to the evolution of the stellar angular momentum. 
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1.1. Continuum-driven Winds from Rotating Stars 

A complete analysis of the dynamics of LBV eruptions is 
quite complicated. In general, one should solve the Navier- 
Stokes equations for the wind, coupled to the radiative trans- 
fer. Furthermore, as we argue that the star is rotating near 
breakup, its oblique shape must be taken into account. Such 
an analysis is amenable only by a numerical treatment. Since 
we wish to understand the underlying physics, this path will 
not be taken here. 

In order to encapsulate the relevant phenomena within a 
simple analytical model, we make the following simplifying 
assumptions. 

I. The star is spherical. 

II. There is no latitudinal motion (vg <» r ), 

Under these approximations, the model is integrable. 

We first note that the typical wind velocity is of the order of 
magnitude of the escape velocity, which is much larger than 
the speed of sound. This implies that the kinetic energy den- 
sity is significantly higher than the internal energy of the gas; 
mechanical pressure and viscosity will therefore be neglected. 

We focus on two fields: the velocity of the wind v and the 
radiative flux F, and begin with conservation of momentum: 



Here £ f = p(g gra , 
per unit mass, with 
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+ v-Vjv = ^]f. (1) 
grad) is the sum of the different forces 
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where g = GM/R 2 and k is the opacity. At a steady state 
the radial component of eq. (1) is 1 
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Conservation of angular momentum implies that 



and hence 
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where 



n = LJy/RZ/GM 



(3) 



(4) 



(5) 
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is the spin in units of the classical breakup limit. 

The radiation field loses energy as it accelerates the wind, 
lifting it over the gravitational pull. Conservation of energy 
then implies 



V • F = pF 



Under these approximations the effective gravity is radial. 



(8) 



which reduces to 



r 2 dr 



r 2 F r 



-pF r v r , 



(9) 



assuming axial symmetry and neglecting latitudinal motion. 
Conservation of mass, 



! + V>v) = 0, 



(10) 



closes the set of equations. At a steady state with the assumed 
symmetries, it reduces to 



For convenience, we define 
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Awr pv r 
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(12) 



where AI is the mass loss rate and is therefore negative. 2 Plug- 
ging this into eq. (9) and integrating gives 



F r = 



L 



Aitr 2 



,-m/(fi)(l-H/r) 



x(0), 



(13) 



where the luminosity L was introduced by dimensional anal- 
ysis in order for the constant of radial integration x(&) to be 
dimensionless, and 

(14) 

AttcR 

is the "photon tiring number". 3 

Using von Ziepel's theorem (von Zeipel 1924) with the sim- 
plification of having just a radial flux, we can determine the 
radiative and mass fluxes at the base of the wind: 
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The requirement J F • dS = L gives 
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(15) 



(16) 



Unlike line driven winds, the mass flux expected in a SED 
continuum driven wind is inherently local. This is because the 
critical surface of the wind, where the gravitational and radia- 
tive forces balance, depends only on the size of the inhomo- 
geneous structure formed by radiative hydrodynamic instabil- 
ities. Since the hydrostatic scale height, over which these 
instabilities operate, is much smaller than the stellar radius, 
mass loss is determined by local conditions. As the local at- 
mospheric structure cannot distinguish between gravity and 
the centrifugal force, the mass loss will depend on the local 
effective gravity. 

Generalizing the local mass flux derived by Shaviv (2001b), 
we have 

W 

* = — (F - F crit ) , (17) 



where 



Fcrit 



'edd 



(1 - n 2 Sil 



(18) 



2 This choice was made to keep / (8) = 1 in the case of no rotation. 

3 This definition is slightly different from the original one defined by 
Owocki & Gayley (1997). Here the mass loss is compared with the Edding- 
ton luminosity, rather than the actual one. 
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This is analogous to the local wind flux in SED accretion 
disks (Dotan & Shaviv 2011). Here F e dd = gc/ 'n is the "Ed- 
dington flux", for which the radiative flux F balances gravity 
in the case of no rotation; F CI i t is the equivalent flux when 
rotation is added. Last, W is the wind constant, which en- 
capsulates mostly geometrical features described by Shaviv 
(2001b), but may have a weak dependence on Y = F j i^dd- 

Owocki et al. (2004) have shown that this dependence is 
stronger if the inhomogeneous structure comprising the un- 
stable atmosphere is power law dominated by smaller scales. 
However, as we shall see below, this merely scales the mass 
loss by a constant factor since Y is latitudinally independent. 
Note that a weak latitudinal dependence of the mass loss 
should arise from the latitudinal dependence of v s at the base 
of the wind. However, v s is proportional to y/T which itself 
is a very weak function of the optical depth. 

Neglecting the aforementioned weak dependences, the flux 
depends on Y as 



<& = pv oc F — F c 



oc 



i + -n 2 ) x (o). 



Under the above approximations, a comparison with eq. (12) 
reveals that 

f(e) = x (0). (20) 



A surface integration of <&, 



M 



gives 



m= -W-SZL 

2 cv s 



Y - 1 



* dS, 



n 2 = w l r - i 



(21) 



(22) 

where W is the scaled wind constant and v 2 sc = 2GM / R. 

An explicit form of F r (eq. 13) allows a direct integration 
of v r (eq. 9): 
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Neglecting the velocity at the base of the wind gives 
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(24) 
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Note that the wind approaches its terminal velocity after 
traversing just a few stellar radii: 



Voo 



> 1 



R 

r 



(26) 



This justifies our second assumption, that vg <C v r . In prin- 
cipal, one could argue that even if the star were spherical, ra- 
diative diffusion would have given rise to latitudinal radiation 



and velocity components. However, since most of the accel- 
eration takes place close to the surface, these components are 
bound to be small. 

1.2. Self Spin-up by Continuum-driven Winds 

As the star blows wind, it loses angular momentum. How- 
ever, specific angular momentum and thus the dimensionless 
spin O may increase, as we now show. 

We start with conservation of angular momentum. The an- 
gular momentum reduction of the star is the angular momen- 
tum taken by the wind: 



where 



and 



2„2 



uiAIR a 



R 2 sin 2 6 * • dS = -luMR 2 



(27) 



(28) 



X sin 3 6d6 



= -u:MR 2 ai 



(29) 



Here a g is the stellar radius of gyration. a w is an effective 
radius of gyration of the wind, defined in the equation above. 
Dividing eq. (27) by l s , we find 



3M 1R 
2M 2R 



CI 



a 2 g M' 



(30) 



In the upper part of the main sequence, where radiation pres- 
sure dominates, one roughly has that R oc VM, such that 
2R/R w M/M. Eq. (30) can now be integrated to give 



% (o? w 7\ M f 



(31) 



This assumes that over the integration interval, a w remains 
constant. 4 The condition for self spin-up is therefore 



(32) 



For a high load wind o? w rj 2/5, and the star spins down. 5 
However, at low mass loss rates, corresponding to m <C 1, 
the flux is insufficient to push the wind over the effective po- 
tential at equatorial latitudes. The wind reaches infinity only 
for angles between the pole and 



2/W-l 



+ G(m 3 / 2 ). 



(33) 



We expect the rest of the wind to stagnate and fall back, thus 
taking no net angular momentum. Therefore, 



a 2 = 



X sin 3 6d6 



(34) 



9 4 

max 



0(6\ 



6 ) 

max ) 



4 This assumption is applicable to both the present day wind with its low 
mass loss rate, and high load winds at a steady state. Eruptions are treated 
differently in the following section. 

5 Massive stars typically have c? g 



■ 0.1 (Motz 1952). 
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One should note that the photon tiring number should be com- 
puted using the mass loss at the surface of the star, not at infin- 
ity. This quantity relates to the observed ejecta by truncating 
the surface integration of <fr (eq. 21) at 8 maK : 



moo =m x sm ' 
lo 



(35) 



2/W-l 



(1 - §f! 2 ) 







Using equations (33) - (36) one finds 

2/VV - 1 



m c 



2f2 2 (l-|Q 2 ) 



(36) 



1.3. Spin Evolution During Eruptions 



In the previous section, we studied the effects of continuum 
driven winds on the evolution of angular momentum assum- 
ing a constant luminosity. As we shall now see, an abrupt 
change in luminosity has significant consequences regarding 
the evolution of the spin of the star. 

A sudden increase in luminosity can be understood as 
a result of an atmospheric phase transition, where the at- 
mosphere becomes porous and the effective opacity drops 
(Shaviv 2001b). Consequently, the radiative flux at the sur- 
face becomes greater than the incoming convective flux. In 
order to compensate for the imbalance, the star must contract 
while radiating the difference in binding energy U : 6 



• I dU ■ 



dR 



(37) 



where L emp and are the luminosities during and before 
the eruption, respectively. A more convenient choice of pa- 
rameters is Lemp and A = 1 — L; n j t /L emp \ henceforth, L e rup 
will be denoted by L. 

In order to keep this discussion as generic as possible, we 
parametrize the binding energy as 



U 



R 



B 



(38) 



and assume a constant B, which is generally the case. The 
time scale for contraction is then 



5 = TkhBA" 
R 



(39) 



where Tkh = GM 2 /LR is the Kelvin-Helmholtz time scale. 
Plugging this into eq. (30) gives 



X (T BV x + (< 3 V 7 - Q 



(40) 



The contraction terminates when the temperature build up at 
the core becomes sufficient to generate the required luminos- 
ity through nuclear reactions. 

In order to estimate how long does the star need to con- 
tract, we assume that it behaves homologously and consider 

6 Because of its low density, the atmosphere alone does not have sufficient 
gravitational binding energy. Note that the increase in luminosity is associ- 
ated only with the release of binding energy due to contraction, which will 
eventually stop once nuclear reactions set in. 



the scaling of the luminosity. If the specific energy production 
can be written as e cx p p T q , we get 

L = J epdV cx pl+PT^R 3 . (41) 

In the following, we assume that p = 1, which is suitable 
for most nuclear reactions, Hydrogen burning included. The 
temperature scales as T cx P/p when gas pressure dominates, 
and T cx P 1 / 4 when radiation pressure dominates. If the dy- 
namical time scale, Td yn = 1/ y/Gp, is much shorter than the 
contraction time, one can safely assume mechanical equilib- 
rium. 7 Hydrostatics then give that P cx M 2 /R 4 . We thus find 
L cx M 2+q/ " R~ 3 ~ q with v = l(or 2), corresponding to the 
case where gas (or radiation) pressure dominates. 

Because q is typically very large, the star needs to contract 
only by a relatively small amount 



AR C , 



R 



1/(3+9) 



AAA 



(2+q/u)/{3+q) 



1 ~~M~) 



(42) 



(43) 



The contraction will therefore take place for a duration 

Jcont — / — •; JKH»A — . 

J R tt 

2. PART II - APPLICATION TO r?-CARINAE 

We now show how the theory developed in §1 applies to 
the particular case of 77-Carinae. In §2.1, it is used to estimate 
the shape of the SED outflow ejected during the Great Erup- 
tion. We argue that the Homunculus' shape is a natural con- 
sequence of a star rotating near breakup with SED continuum 
driven winds. In §2.2, we show that the present day wind is 
sufficiently concentrated towards the poles to spin up the star 
towards critical rotation, if it is the dominant mode of mass 
loss. In §2.3, we show that also the formation of an equatorial 
skirt is expected given the stellar parameters of ^-Carinae. 

We consider nominal values of M = 12OM , R = 15Oi?0 
for the mass and radius. Since the age of the star is unknown, 
we use an intermediate chemical composition corresponding 
to k = 0.3 cm 2 g _1 for the opacity. As advocated, the star is 
assumed to be near critical rotation. 

2.1. The Homunculus as a Continuum-driven Wind from a 
Rotating Star 

We assume that during the Great Eruption, the sus- 
tained mass loss rate had been of the order of M = 
— 0.1M Q year _1 . 8 Given the aforementioned nominal stel- 
lar parameters, the corresponding photon tiring number is 
m w 0.5. Assuming critical rotation, the size of the Ho- 
munculus is reproduced by a value of T w 1.6, which in turn 
implies a scaled wind constant of W ~ 0.4. This value of V 
corresponds to a luminosity L = rL e dd ~ 10 7 L Q at the base 
of the wind. 

At infinity, the modeled luminosity must account for photon 
tiring: 



F • dS = L 



47T 



-mx(0). 



e -^~'x{0) ~ 0.5L. (44) 



For comparison, the observed luminosity at the time was in 
the range of L obs « 1 - 3 x 10 7 i© (Smith & Frew 201 1). 

7 For the case of r)-Carinae, T,j yn f» 2 weeks. 

8 More about this choice will be discussed below. 
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Figure 1. The latitudinal dependence of the radial velocity of the 
wind. The dots correspond to the values observed by Smith (2006). 
The solid line represents the model described in §1.1 for the nominal 
stellar parameters and T « 1.6 at critical rotation. The dashed line 
corresponds to the same model with f2 R2 0.9 - the velocity profile 
expected after a few years of mass loss at a rate of M = — O.lMQyear - 1 . 



Figure 1 shows the latitudinal dependence of the wind ve- 
locity predicted by eq. (24) versus observed values. 

2.2. Spin-up from the Present-day Wind 

The present day mass loss rate is of the order of 
lO -3 M0year -1 (Davidson & Humphreys 1997), which im- 
plies ps 3 x 10~ 2 (eq. 36). In this mode of mass loss, it 
would take Tj-Carinae about 2 x 10 3 years to reach again crit- 
ical rotation, for the above nominal parameters. The maximal 
mass loss rates able to satisfy eq. (32) as a function of spin 
can be read off Figure 2. 



2.3. The Formation of the Equatorial Skirt 

Any attempt to model the dynamics behind the Homunculus 
nebula cannot be complete without explaining the formation 
of the equatorial skirt. As it turns out, the formation of such a 
skirt is a generic property of SED eruptions at critical rotation. 

Consider the equation of conservation of angular momen- 
tum for the contracting star (eq. 40). If the star erupts while 
critically rotating and 



M 
2kh 



> \M\ 



(45) 



the equator becomes unbound, and a skirt is formed at a rate: 



.9 



Ma 



A. 



-Br 



M 



cvr, 



M 



(46) 



Note that the mass loss to the skirt is neglected with respect to 
the mass loss to the wind; this is justified by self consistency. 

For a critically rotating n = 3 polytrope, which is a reason- 
able description of ?y-Carinae, one has that 



B 



1 2 

2°" 



(47) 



where (3 = P gds /Pt t (see, for example, Clayton 1968). 

Shaviv (2001a) has shown that the radiative hydrodynamic 
instabilities, which are responsible for the atmospheric phase 

9 This can be seen by replacing Q/Q by M^/M or in eq. (40). 



Figure 2. The spin-up condition (eq. 32) is satisfied for values of spin and 
mass loss rate below the line. 




0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 

P 

Figure 3. The mass of the equatorial skirt, as a function of M and ft = 
P, gas /Pi tot- The values given here should be considered as lower limits, as the 
relevant radius of gyration of the outer shell participating in the dynamics is 
much larger than that of the whole star, see §B. 

transitions, are excited for 0.15 < f3 < 0.5, depending on 
different conditions. In this range of j3, one has v ps 2. 

Comparing the luminosities computed in the previous two 
sections, we see that during the eruption, the luminosity in- 
creased roughly by a factor of 5, and so A ps 4/5. 

The last piece of information required to understand the 
contraction is the factor q, relating the specific energy produc- 
tion and temperature by e oc T q . Since the temperature at the 
core of 77-Carinae is expected to be of the order of 5 x 10 7 K, 
the energy production is due to CNO burning, which at this 
temperature corresponds to q ~ 15 (e.g., Bethe 1939). 

We can now estimate (eq. 43) the contraction time to be 3 
to 10 years, depending on the value of j3, and so the skirt's 
mass is of the order of 0.1 Mq. Figure 3 shows the predicted 
mass of the skirt for different values of j3 and M. 

Throughout this analysis, we have implicitly assumed that 
the whole star is instantly affected by the loss of angular mo- 
mentum to the wind, and thus the corresponding a 2 g is con- 
stant. However, during time periods shorter than the time 
scale of angular momentum diffusion, only an outer shell 
plays a role in the dynamics. In §B, we show that the dif- 
fusion time scale is indeed comparable if not longer than the 
contraction time. As a result, the effective radius of gyration 
is larger, and so the mass of the skirt computed from eq. (46) 
is in fact a lower limit. We conclude that during the initial 
stages of contraction, a skirt is likely to have formed even for 
relatively high mass loss rates. 
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3. DISCUSSION AND SUMMARY 

The twin lobe structure of the Homunculus nebula of 77- 
Carinae strongly suggests that either fast rotation or binary 
interaction played an important role in the process of mass 
loss. In this work, we have chosen to pursue the idea that it is 
near critical rotation that sculpted the Homunculus, and that 
the mass loss was in the form of a continuum driven wind. 
The formation of the equatorial skirt is then associated with 
the evolution of angular momentum before the system settled 
in a steady state. 

As we have seen in §1.1, the structure of such winds can 
be obtained analytically, if one allows several simplifying as- 
sumptions. In particular, it was assumed that there is no latitu- 
dinal interaction between mass and radiative flux elements on 
different radial trajectories. This assumption later proved to 
be self consistent, as it turns out that most of the acceleration 
of the wind takes place near the surface. 10 Namely, the rele- 
vant radial scale is smaller than the latitudinal one; this acts 
to suppress lateral fluxes. 

By further simplifying the geometry and assuming a spher- 
ical star, the set of equations describing all conserved quan- 
tities became both algebraically closed and integrable. Ob- 
viously, this assumption cannot hold since rotation breaks 
spherical symmetry. Even so, we find that this simple model 
captures the gross features of the system, as depicted in fig- 
ure 1. Mainly, a twin-lobed structure emerges as a generic 
property of winds from such stars. The main ingredient re- 
quired in order to form this structure is that of gravity darken- 
ing, as was already pointed out by Owocki et al. (1998) and 
Maeder (1999) for line driven winds. Nevertheless, we do ex- 
pect the oblate geometry to introduce distortion. For example, 
the predicted velocity may be lower near the equator, because 
the wind is launched from larger radii. This may in fact im- 
prove the agreement with the observations. 

For the specific case of 77-Carinae (§2.1), our model pro- 
vides a rough description of the shape and size of the Ho- 
munculus nebula. In particular, continuum driven winds natu- 
rally explain the large wind momentum to photon luminosity 
observed in the Homunculus. Such high ratios cannot be ex- 
plained by line driven winds. 

Note however, that within the approximations and the nom- 
inal parameters used, the typical mass of the Homunculus 
obtained is somewhat lower than the value recent observa- 
tions imply (Smith et al. 2003), and the predicted luminosity 
is lower by a factor of ~ 3. These discrepancies may be re- 
duced by a full numerical analysis. 

As a side note, we show in § A that the angular momentum 
gain due to a general spin-orbit coupling is correlated with the 
change in semi-major axis. As a consequence, in light of the 
high eccentricity of the orbit and stellar parameters assumed, 
the angular momentum 77-Carinae may have received from in- 
teractions with its companion were insufficient to drive it to 
critical rotation. 

The fact that the rotation is very close to critical may seem 
like fine tuning, but it is quite possible that this is the nat- 
ural state of the system. Since the star cannot sustain high 
load winds for long periods of time, it is reasonable that a low 
mass loss rate, as ?7-Carinae currently exhibits, is in fact the 
dominant mode of mass loss. Motivated by this reasoning, we 

10 Eq. (24) implies that the winds of r/-Carinae reach half their terminal 
velocity already at r ~ 1.1 5R. For comparison, line driven winds reach half 
their terminal velocity further out, at r m 1.73-R (e.g., Lamers & Cassinelli 
1999). 



analyzed in §1.2 the explicit solutions for the wind velocity 
at low mass loss rates. We found that at equatorial latitudes, 
the wind stagnates and falls back to the surface, given that the 
star is rotating fast enough. The mass loss from pole centered 
winds, which do reach infinity, dominates over the loss of an- 
gular momentum and so the star spins up. 1 1 This is another 
aspect of continuum driven winds which is absent from the 
line driven winds model. 

We predict (§2.2) that in about 2 x 10 3 years at the current 
mass loss rate, 77-Carinae will be back at critical rotation. This 
implies that if the previous large eruption cycle and subse- 
quent wind were similar to the Great Eruption and the present 
wind, then the previous large eruption of 77-Carinae must have 
taken place at least a few thousand years ago. 

We then showed in §1.3 that at the onset of LBV eruptions, 
the contraction that follows the atmospheric phase transition 
is fast enough to render equatorial latitudes unbound, if the 
star is critically rotating. The formation of an equatorial skirt 
is seen to be a generic consequence of the model as well. 
For the stellar parameters of jy-Carinae (§2.3), we predict the 
skirt's mass to be of the order of 0.1 Mq. 

The effect of angular momentum diffusion, which is dis- 
cussed in §B, complicates the estimation of the mass of the 
skirt. Taking it into account requires a more accurate knowl- 
edge of the mass loss to the wind, and specifically its time 
dependence. Essentially, the question is what happened first 
- did nuclear reactions kick in to stop the contraction or did 
the loss of angular momentum manage to diffuse inwards and 
render the formation of a skirt unnecessary? We await the 
appearance of new observational data regarding the time de- 
pendence of the Great Eruption. 
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APPENDIX 

A. A BOUND ON THE ANGULAR MOMENTUM GAINED FROM 
SPIN-ORBIT COUPLING 

The above analysis assumed that there is no additional ex- 
change of angular momentum. However, 77-Carinae is part of 
a binary system, and one might suspect that spin-orbit cou- 
pling plays a dominant role in its spin evolution. Given the 
orbital parameters, we argue that such interactions do not of- 
fer a significant contribution to the spin, and can therefore be 
neglected. 

The basic idea is that the change in the spin is correlated 
with a change in the semi-major axis; any angular momentum 
conserving interaction cannot change one without effecting 
the other. We show that even if the original binary orbit had an 
extremely large semi-major axis, the evolution to the present 
orbit could not have spun up 77-Carinae significantly. 

Let us then consider a general force f representing the in- 
teraction. Assuming the interaction is important only at short 
distances, 12 it may be approximated as an impulse during pe- 
riastron passages. 

1 1 By "spin-up" we refer to an increase of Q = ui yj B? /GM. 

12 For example, tidal interactions give that f oc r~ 7 (Hut 1981). 
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Due to conservation of total angular momentum, the compar- 
ison between the spin angular momentum, £ s , and the orbital 
angular momentum, £ a , gives 



£ = — [for + rr/j.9 



(Al) 



where r(8) describes the orbital motion. We further neglect 
the term rr/i# since r — > at periastron. In order to reduce 
the clutter of notation, £ will be henceforth used to denote £ s . 
During one orbit, the change in the angular momentum is 



Si: 



£dt 



dt. 



(A2) 



where the subscript p denotes values calculated at perias- 
tron. 13 The change in the semi-major axis is due to the loss of 
orbital energy, E = —GMmjla, hence 
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On the other hand, the change in energy is 



SE= / f • dr 



f rdQ w -£J p / dt 



(A3) 
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One should note that the periastron separation, r p = a(l — e), 
is approximately constant during orbital evolution at high ec- 
centricities. Since the interaction takes place at periastron, the 
evolution of the orbit is such that the periastron remains fixed 
while the apastron, r a , progressively shrinks. Although the 
orbit may change during periastron passages when the force 
is applied, the periastron itself must remain part of it. This 
implies that the angular velocity at periastron passages, 



r p 3 G(M + m)(l + e), 



(A5) 



changes very slowly and is assumed to be constant as well. 

On a long time scale, the ratio between time derivatives of 
the semi-major axis and the spin is the ratio between the cor- 
responding incremental changes every periastron passage: 
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This expression can be integrated to give 

t 
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Clearly, the maximum possible change in the angular momen- 
tum is 

GMm 



A£ n 



29 p a 



(A8) 



Namely, the upper bound depends only on the current value 
of the semi-major axis. 

Another way of looking at the orbital evolution is global. 
We can compare the angular momentum of two systems hav- 
ing the same periastron, r p , but one having r a — s- oo, and the 
other an r a of the present orbit. In contrast to the previous 
derivation, a global view cannot explain why the perisatron 
distance remains approximately constant during the evolution. 

13 We use 8 to denote changes over one orbit, and we will use A for 
changes over the whole evolution. 



The general expression for the angular momentum of a Ke- 
plerian orbit is 



G 



M 2 m 2 
M + m 



r p (l + e). 



(A9) 



We define A£ omax = £(e) — £(e —> 1) and expand the angular 
momentum for small values of 1 — e. We then obtain, to the 
leading order, 

|A£ , max | < (A10) 
26 p a 

in agreement with eq. (A8). 

This change of angular momentum should be compared 
with the present value of the angular momentum of the star, 
while assuming it is rotating near breakup. For the present 
nominal stellar parameters, including a period of 5.5 years 
and e = 0.9, we find that 



A£ m J£ s ~ 5%, 



(All) 



and thus interaction with the companion star could not be held 
responsible for the high spin of 77-Carinae. 

B. ANGULAR MOMENTUM DIFFUSION 

The last theoretical calculation required for the understand- 
ing of angular momentum evolution in mass losing stars is an 
analysis of angular momentum diffusion. Until now, we have 
implicitly assumed that the whole star immediately reacts to 
the mass loss. However, the reaction time cannot be shorter 
than the diffusion time scale Tdiff. 

When considering processes lasting for time scales compa- 
rable to Tdiff, the star cannot be thought of as a rigidly rotating 
body. Since only an outer shell participates in the dynamics, 
the effective radius of gyration is much larger than that of the 
entire star. 

The goal of this section is to obtain a lower bound for Tdiff, 
the diffusion time scale associated with the homogenization of 
such a shell. As we shall see, the angular momentum diffusion 
time in massive stars ranges from a fraction of year for the 
outer parts of the envelope, to several years to homogenize 
the whole star. 

Since Tdiff ~ R 2 /D, we seek an upper bound on the dif- 
fusion coefficient D. We note that the most efficient diffu- 
sion mechanism in stars is that of convective mixing (Taylor 
1915). 14 We also note that convectionis always excited above 
some critical, sub-Eddington luminosity (Jossefa/. 1973), 
implying that stars with continuum driven winds are neces- 
sarily convective. This motivates the application of mixing 
length theory (MLT, Vitense 1953; Bohm-Vitense 1958) to 
angular momentum transport in a SED state. 

Under MLT, the diffusion coefficient is 



V 



hv c 
~3~ 



(Bl) 



where v c is the convective velocity and h is the pressure scale 
height, playing the role of the mean free path. 
MLT also gives that the convective velocity is 



L 



\~4tt 
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(B2) 



14 Other mechanisms such as shear induced turbulence or meridional cir- 
culation are important when convection is absent (e.g., Zahn 1992). 
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Figure 4. The solid line depicts the squared radius of gyration of a shell 
located between a certain distance from the center, r, and the stellar surface, 
R. The dashed line plots the dimensionless function T describing the typical 
time it takes to convectively diffuse angular momentum throughout the shell. 
These functions are plotted for an n = 3 polytrope. 

where a ~ c p p/v^ «3/2 and c p is the heat capacity per unit 
mass. Since v c is obviously limited by the speed of sound, 
eq. (B2) should be regarded as an upper bound. 

In order to compute the scale height, we model the star as 
an n = 3 polytrope 15 (Eddington 1931): 



h = 



dr 



= R 



1 
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The diffusion coefficient is then 



V = 
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and the typical diffusion time between a radius r and the sur- 
face of the star is approximately given by 




(B7) 



For the stellar parameters of ?7-Carinae, the diffusion time is 

T m ~ 10 years x (B8) 
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15 For completeness, we briefly review the polytropic approximation: the 
pressure and density are conjectured to be related by p = p c <t> n and P = 
P c <fi n + for some value of n. The hydrostatic equation, VP = — pV$, then 
reads in + 1)P C = — p c V<E>, and Poisson's equation for the gravitational 
potential, V 2 $ = 4-trGp, becomes the Lane-Emden equation (Lane 1870; 
Emden 1907): 



1 d 



0", (B3) 
"^TrGp^n + lJ-ip- 1 . 



where spherical symmetry is assumed and £ 
£l is the value for which = 0; for n = 3, a numerical computation 
gives £i as 6.9 . Note that £i is related to the physical radius through 



The radius of gyration of the outer shell is given by 



Hi 2 // 1 £ 2 « 



(B9) 



The functions T and a 2 h are plotted in figure 4. 
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